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PROCEDE D'ESTIMATION DE DISTANCE POUR UN MOBILE SOUMIS A 
DES CONTRAINTES DYNAMIQUES DE PARCOURS 

L'invention concerne la navigation de terrain pour un mobile 

5 soumis ^ des contraintes de parcours variant au cours du temps, tel qu'un 
aeronef limite en taux de montee, la limite pouvant §tre n6gative, et 6voluant 
au-dessus d'une zone de terrain prSsentant des reliefs ou des obstacles 
menagants proches ou superieurs d son altitude de vol. 

Divers systfemes ont 6te d^veloppes pour pr6venir T^quipage d'un 

10 a6ronef d'un risque de collision avec le sol. Certains, tels que les systfemes 
TAWS (acronyme de Texpression anglo-saxonne : Terrain Awareness and 
Warning System"), font une pr6vision de trajectoire a court terme pour 
ra6ronef ^ partir des informations de vol (position, cap, orientation et 
amplitude du vecteur Vitesse) fournies par les equipements du bord, la 

16 placent en situation par rapport a une carte de la r§gion survol6e extraite 
d'une base de donnees d'6levation du terrain accessible du bord et emettent 
des alarmes a destination de I'equipage de I'aeronef a chaque fois que la 
trajectoire previsible § court terme entre en collision avec le sol. Ces 
systemes TAWS agrementent leurs alarmes de recommandations 

20 rudimentaires du genre "Terrain Ahead, Pull up". Certains d'entre eux 
donnent 6galement des informations sur le niveau de risque de collision que 
font encourir les reliefs et les obstacles environnant I'aferonef sous la forme 
d'une carte pr6sentant les reliefs ou les obstacles du terrain survole en 
strates de couleurs diff^rentes. Cependant, cette carte de risques de collision 

25 avec I'environnement tient uniquement compte des altitudes du relief 
relativement d la position du mobile et ne prend pas en compte de Texistence 
ou non d'une trajectoire r6aliste permettant de rejoindre les zones affich6es. 

Pour satisfaire ce besoin de connaTtre les points du terrain survole 
restant accessibles aprds une manoeuvre d'6vitement d'un relief ou d'un 

30 obstacle au sol, la carte de risque de collision avec I'environnement doit 
afficher uniquement les zones pour lesquelles 11 existe un chemin possible 
depuis la position courante du mobile. La realisation d'un tel affichage passe 
par Tassociation d'une metrique d une carte du relief tiree d'une base de 
donnees d'elevatipn du terrain. 

35 Une m6thode connue pour associer une metrique a une carte du 

relief tir6e d'une base de donnees d'61evation du terrain d maillage r6gulier 
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de la surface terrestre ou d'une partie de celle-ci, consiste ^ consid^rer la 
carte pr6sentant le relief d partir de valeurs d'altltude figurant, avec les 
coordonn6es g§ographiques : latitude et longitude des points de mesure, 
dans les elements de la base de donn^es delevation du terrain comme une 

5 image dont les pixels sont les valeurs d'altitude des points de la base de 
donnSes d'el6vatton du terrain illustrSes dans la carte avec, pour 
coordonnSes abscisse et ordonnSe au sein de I'image, les coordonnees 
geographiques latitude et longitude de ces points figurant dans les elements 
de la base de donn6es d'^levation du terrain et a faire appel a une 

10 transform6e de distance operant par propagation pour estimer des distances 
au sein de cette image. 

Les transformees de distance operant par propagation 6galement 
connues sous Tappellation "transformees distance de cfianfrein" ("chamfer 
distance transform" ou "chamfer euclidean distance transform" en anglo- 

15 saxon) deduisent la distance d'un pixel dit pixel but par rapport S un autre 
pixel dit pixel source, des distances prec^demment estimSes pour les pixels 
de son voisinage, par un balayage des pixels de I'image. Le balayage permet 
d'estimer la distance d'un nouveau pixel but par rapport au pixel source par 
recherche du trajet de longueur minimale aliant du nouveau pixel but au pixel 

20 source en passant par un pixel interm^diaire de son voisinage dont la 
distance a d6jd 6t6 estim^e, la distance du nouveau pixel but d un pixel 
interm6diaire de son voisinage dont la distance a deja 6te estim6e 6tant 
donn^e par application d'un masque de voisinage commun^ment appeld 
masque de chanfrein. 

25 Une transform6e de distance de ce genre a 6t6 proposee en 1986 

par Gunilla Borgefors pour estimer des distances entre objets dans une 
image num6rique, dans un article intitule :" Distance Transformation in Digital 
Images." et paru dans la revue "Computer Vision, Graphics and Image 
Processing", Vol. 34 pp. 344-378. Un des interets de ces transformees de 

30 distance par propagation est de r6duire la complexite des calculs d'une 
estimation de distance en autorisant I'emploi de nombres entiers. 

Pour selectionner le trajet de longueur minimale donnant 
I'estimation distance, une transformee de distance par propagation doit tester 
tous les trajets possibles. Cette obligation se traduit par une contrainte de 

35 r6gularit§ impos^e d I'ordre de balayage des pixels d'une image. G. 
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Borgefors propose, pour satisfaire cette contrainte de regularit6, de balayer 
les pixels d'une image deux fois cons^cutivement, dans deux ordres inverses 
Tun de Tautre, qui sent soit Tordre lexicographique, Timage etant analysee de 
gauche d droite ligne par ligne et de haut en bas, et I'ordre lexicographique 

5 inverse, soit Tordre lexicographique transpose, rimage ayant subi une 
rotation de 90**, et Tordre lexicographique transpos6 inverse. Elle propose 
Sgalement d*adopter un masque de chanfrein de dimensions 3x3 avec deux 
valeurs (3, 4) de distances de voisinage ou de dimensions 5x5 avec trois 
valours (5, 7, 1 1) de distances de voisinage. 

10 Les transform6es de distance operant par propagation sont d^ja 

employ6es dans le domaine de la navigation terrain pour robot. Dans ce 
cadre, il est connu d'utiliser la transformee de distance de G. Borgefors avec 
une contrainte statique consistent a attribuer, de maniere autoritaire, une 
distance infinie ^ un point en analyse lorsqu'il apparaTt qu'il fait partie de 

15 reliefs ou d'obstacles S contourner repertories dans une m6moire de zones 
interdites de franchissement. cela de manifere d 6liminer, de I'ensemble des 
trajets testes lors d'une estimation de distance, ceux passant par les reliefs 
ou obstacles que le robot doit contourner. Cependant, une transformee de 
distance operant par propagation utilisde avec une contrainte statique dans 

20 le cadre de la navigation terrain pour robot, ne convient pas d la navigation 
terrain pour a§ronef pour lequel la menace pr^sentSe par un relief ou un 
obstacle au sol depend du profil vertical de sa trajectoire. 

La pr^sente invention a pour but un precede d'estimation des 
25 distances des points d'une carte extraite d'une base de donnees d'elevation 
du terrain par rapport a un point de reference employant une transformee de 
distance operant par propagation, avec une contrainte dynamique Evolutive 
au cours du temps convenant a la navigation terrain pour un aeronef ayant 
une trajectoire a profil vertical impose. 

30 

Eile a pour objet un precede d'estimation des distances des points 
d'une carte extraite d'une base de donnees d'6l6vation du terrain, pour un 
mobile soumis d des contraintes dynamiques lui interdisant certaines zones 
de la carte dites zones interdites de passage dont la configuration varie en 
35 fonction du temps de parcours du mobile. La base de donn6es d'§l6vation du 
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terrain renferme un ensemble de points reper^s par une altitude, une latitude 
et une longitude maiilant le terrain d'^volution du mobile. Le proc6d§ met en 
oeuvre une transfonm6e de distance op6rant par propagation sur Timage 
constitu6e par les 6lements de la base de donnees d'el6vation. du terrain 

5 cx)rrespondant a la carte, dispos6s en lignes et colonnes par ordres de 
valeurs de longitude et de latitude. Cette transformee de distance estime les 
distances des differents points de I'image par rapport a un point source place 
d proximity du mobile en appliquant. par balayage, un masque de chanfrein 
aux differents points de Timage. L'estimation de distance d'un point, par 

10 application du masque de chanfrein a ce point dit point but s'effectue en 
rfepertoriant les diff6rents trajets allant du point but au point source et passant 
par des points du voisinage du point but qui sont recouverts par le masque 
de chanfrein et dont les distances au point source ont 6t6 pr^atablement 
estim^es au cours du meme balayage, en determinant les longueurs des 

15 differents trajets repertories par sommation de la distance affectee au point 
de passage du voisinage et de sa distance au point but extraite du masque 
de chanfrein, en recherchant le trajet le plus court parmi les trajets 
repertories et en adoptant sa longueur comme estimation de la distance du 
point but. Initialement, en debut de balayage, une valeur de distance 

20 superieure d la plus grande distance mesurable sur Timage est attribuee d 
tous les points de Timage sauf pour le point source, origine des mesures de 
distance, auquel est affectee une valeur de distance nulle. Le precede est 
remarquable en ce que les longueurs des trajets repertories, lors de 
rapplication du masque de chanfrein d un point but, en vue de la recherche 

25 du trajet le plus court, sont traduites en temps de parcours pour le mobile et 
en ce que les trajets repertories dont les temps de parcours pour le mobile 
sont tels que le point but appartiendrait a une zone interdite de passage au 
moment ou le mobile I'atteindrait, sont exclus de la recherche du trajet le plus 
court. 

30 

Avantageusement, lorsque le mobile est un aeronef ayant un profil 
vertical de vol e respecter determinant revolution de son altitude instantanee, 
on associe, aux longueurs des trajets repertoriees, les valeurs previsibles 
des altitudes Instantanees qu'aurait Taeronef en atteignant le point but par 
35 ces trajets tout en respectant le profil vertical de vol impose, et on etimine de 
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ia recherche du trajet le plus court, les trajets r6pertorl§s associ6s d des 
valeurs pr6visibles d'altitude inf6rieures ou egales a celle du point but 
donn6e par la base de donn6es d'6l6vation du terrain et augmentee d'une 
marge de protection. 

5 

Avantageusement, lorsque le mobile est un a6ronef ayant un profil 
vertical de vol impose, I'estimation de distance operee par propagation sur 
rimage constitute des Elements de la base de donnees d'6levation du terrain 
correspondant ^ la carte est doublee d'une estimation de I'altitude previsible 
10 pour I'aeronef au droit des diff6rents points de Timage en supposant qu'il 
suive le trajet le plus court selectionne pour Testimation de distance et qu'il 
respecte le profil vertical de vol impose. 

Avantageusement. lorsque le mobile est un aeronef a profil 
15 vertical de vol impos6 et que Testimation de distance est doublSe par une 
estimation de I'altitude previsible de Taeronef. les altitudes des diff6rents 
points de la carte sont soustraites des estimations des altitudes pr6visibles 
pour I'atronef en ces points pour obtenir des hearts par rapport au sol. 

20 Avantageusement, lorsque le mobile est un aeronef S profil 

vertical de vol impost et que Testimation de distance est doubl6e par une 
estimation de {'altitude prtvisible de I'atronef, les altitudes des differents 
points de la carte sont soustraites des estimations des altitudes prtvisibles 
pour I'aeronef en ces points pour obtenir des 6carts par rapport au sol 

25 affiches sur la carte en strates de couleurs. 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de Timage constitute des elements de la base de donntes 
d'eltvation du terrain correspondant a la carte, en plusieurs passes 
30 successives selon des ordres differents. 

Avantageusement, la transformte de distance par propagation 
balaie les pixels de I'image constitute des titments de la base de donntes 
d'tltvation du terrain appartenant d la carte, en plusieurs passes 
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successives selon des ordres diff6rents et de maniere r6petee jusqu'^i ce que 
les estimations de distance obtenues se stabilisent. 

Avantageusement, la transfomi6e de distance par propagation 
5 balaie les pixels de Timage constituee des 6l6ments de la base de donnees 
d'6l6vation du terrain con^espondant a la carte, en plusieurs passes 
successives selon des ordres differents dont I'ordre lexicographique, Tordre 
lexicographique Inverse, Tordre lexicographique transpose et Tordre 
lexicographique transpose inverse. 

10 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de Timage constitute des elements de la base de donn6es 
d'6levation terrain correspondant a la carte, en une serie de quatre passes 
repetee jusqu'a stabilisation des estimations de distances : 
15 " une premiere passe effectu6e ligne par ligne de haut en bas de 

rimage, cheque ligne etant parcourue de gauche a droite, 

- une deuxi^me passe effectuSe ligne par ligne de bas en haut 
de rimage, cheque ligne ttant parcourue de droite d gauche, 

- une troisitme passe effectute colonne par colonne de gauche 
20 d droite de Timage, cheque colonne 6tant parcourue de haut 

en bas, et 

- une quatridme passe effectu§e colonne par colonne de droite a 
gauche de I'image, cheque colonne etant parcourue de bas en 
haut. 

25 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de Timage constituee des elements de la base de donnees 
d'elevation du terrain correspondant a la carte, en une serie de huit passes 
repetee jusqu'd stabilisation des estimations de distances. : 
30 " une premiere passe effectu6e ligne par ligne de haut en bas de 

I'image, chaque ligne etant parcourue de gauche d droite, 

- une deuxidme passe effectute ligne par ligne de bas en haut 
de rimage, chaque ligne 6tant parcourue de droite d gauche, 



wo 2005/031262 



PCT/EP2004/052078 



7 



une troisi^me passe effectuee colonne par colonne de gauche 
d droite de I'image, chaque colonne 6tant parcoume de haut 
en bas, 

une quatridme passe effectuee colonne par colonne de droite d 
gauche de rimage, chaque colonne etant parcourue de bas en 
haut, 

une cinqui^me passe effectuee ligne par ligne de haut en bas 
de rimage, chaque ligne etant parcourue de droite S gauche, 
une sixi^me passe effectuee ligne par ligne de bas en haut de 
rimage, chaque ligne 6tant parcourue de gauche a droite. 
une septieme passe effectuee colonne par colonne de droite ^ 
gauche de rimage, chaque colonne 6tant parcourue de haut en 
bas, et 

une huitieme passe effectuee colonne par colonne de gauche 
a droite de Timage, chaque colonne 6tant parcourue de bas en 
haut. 

D'autres caracteristlques et avantages de Finvention ressortiront 
de la description ci-apres d'un mode de realisation donn6 d titre d'exemple. 
20 Cette description sera faite en regard du dessin dans lequel : 

- une figure 1 reprfesente un exemple de masque de chanfrein, 

- des figures 2a et 2b montrent les cellules du masque de 
chanfrein illustr6 a la figure 1, qui sont utilis6es dans une 
passe de balayage selon I'ordre lexicographique et dans une 

25 passe de balayage selon I'ordre lexicographique inverse, 

- une figure 3 est un diagramme illustrant les principales etapes 
d'un precede, conforme a Tinvention, pour estimer la distance 
d'un point en tenant compte d'une contrainte dynamique au 
cours de {'application d'un masque de chanfrein, 

30 - une figure 4 est un diagramme illustrant une variante du 

precede d'estimation de la distance d'un point montre a la 
figure 3. et 

- une figure 5 est un diagramme des principales Stapes d'un 
proc6de, conforme d I'invention, d'estimation, par propagation, 

35 des distances de {'ensemble des points d'une carte tenant 
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compte d'une contrainte dynamique et mettant en oeuvre un 
precede d'estimation de la distance d'un point tel que ceux 
montres aux figures 3 et 4. 

5 La distance d'entre deux points d*une surface est la longueur 

minimale de tous les parcours possibles sur la surface partant de Tun des 
points et aboutissant d Tautre. Dans une image formee de pixels repartis 
selon un maillage regulier de lignes, colonnes et diagonales, une 
transformee de distance par propagation estime la distance d'un pixel dit 

10 pixel "but" par rapport a un pixel dit pixel "source" en construisant 
progressivement, en partant du pixel source, !e plus court trajet possible 
suivant le maillage des pixels et aboutissant au pixel but et en s'aidant des 
distances trouvees pour les pixels de I'image d6j^ analyses et d'un tableau 
dit masque de chanfrein r^pertoriant les valeurs des distances entre un pixel 

15 et ses proches voisins. 

Comme montre a la figure 1 , un masque de chanfrein se presente 
sous la forme d'un tableau avec une disposition de cases reproduisant le 
motif d'un pixel entour6 de ses proches voisins. Au centre du motif, une case 
affectee de la valeur 0 repSre le pixel pris pour origine des distances 

20 repertori^es dans le tableau. Autour de cette case centrale, s'agglom^rent 
des cases periph6riques remplies de valeurs de distance non nulles et 
reprenant la disposition des pixels du voisinage d'un pixel suppos6 occuper 
la case centrale. La valeur de distance figurant dans une case p^ripherique 
est celle de la distance separant un pixel occupant la position de la case 

25 peripherique concemee, d'un pixel occupant la position de la case centrale. 
On remarque que les valeurs de distance se repartissent en cercles 
concentriques. Un premier cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places soit sur la ligne 
soit sur ta colonne du pixel de la case centrale sont affectees d'une valeur de 

30 distance D1. Un deuxieme cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places en dehors de la 
ligne et de la colonne du pixel de la case centrale sont affectees d'une valeur 
de distance D2. Un troisi^me cercle de huit cases conrespondant aux huit 
pixels les plus proches du pixel de la case centrale places en dehors de la 
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ligne, de la colonne et des diagonales du pixel de la case centrale sont 
affectees d'une valeur D3. 

Le masque de chanfrein peut couvrir un voisinage plus ou moins 
§tendu du pixel de la case centrale en rSpertoriant les valeurs des distances 

5 d'un nombre plus ou moins important de cercles concentriques de pixels du 
voisinage. II peut 3tre reduit aux deux premiers cercles formes par les pixels 
du voisinage d'un pixel occupant la case centrale ou §tre ^tendu au-deld des 
trois premiers cercles formes par les pixels du voisinage du pixel de la case 
centrale mais il est habituel de s*arreter ^ trois premiers cercles comme celui 

10 represents ^ la figure 1. Les valeurs des distances D1, D2, D3 qui 
correspondent a des distances euclidiennes sont exprimees dans une 
6chelle autorisant Temploi de nombres entiers au prix d'une certaine 
approximation. C'est ainsi que G. Borgefors donne a la distance d1 
correspondant a un 6chelon en abscisse x ou en ordonnee y la valeur 5, a la 

15 distance d2 correspondant a la racine de la somme des carr6s des echelons 
en abscisse et ordonnee ^Jx^ la valeur 7 qui est une approximation de 
5y/2, et a la distance d3 la valeur 1 1 qui est une approximation de sVs . 

La construction progressive du plus court trajet possible allant d 
un pixel but, en partant d*un pixel source et en suivant le maillage des pixels 

20 se fait par un balayage regulier des pixels de Timage au moyen du masque 
de chanfrein. Initialement, les pixels de Timage se voient affecter une valeur 
de distance infinie, en fait un nombre sufTisamment Sieve pour dSpasser 
toutes les valeurs des distances mesurables dans Timage, d Texception du 
pixel source qui se voit affecter une valeur de distance nulle. Puis les valeurs 

25 initiales de distance affectSes aux points but sont mises d jour au cours du 
balayage de Timage par le masque de chanfrein, une mise d jour consistant 
d remplacer une valeur de distance attribute d un point but, par une nouvelle 
valeur moindre resultant d'une estimation de distance faite d Toccasion d'une 
nouvelle application du masque de chanfrein au point but considers. 

30 Une estimation de distance par application du masque de 

chanfrein a un pixel but consiste a repertorier tous les trajets allant de ce 
pixel but au pixel source et passant par un pixel du voisinage du pixel but 
dont la distance a dej^ 6te estimee au cours du meme balayage, ^ 
rechercher parmi les trajets repertories, le ou les trajets les plus courts et ^ 

35 adopter ta longueur du ou des trajets les plus courts comme estimation de 
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distance. Cela se fait en plagant le pixel but dont on veut estimer la distance 
dans la case centrale du masque de chanfrein, en s^lectionnant les cases 
p6riph6riques du masque de chanfrein correspondant d des pixels du 
voisinage dont la distance vient d*§tre mise a jour, en calculant les longueurs 

5 des trajets les plus courts reliant le pixel a mettre a jour au pixel source en 
passant par un des pixels s6lectionnes du voisinage, par addition de la 
valeur de distance affect6e au pixel du voisinage concern6 et de la valeur de 
distance donn6e par le masque de chanfrein, et a adopter, comme 
estimation de distance, le minimum des valeurs de longueur de trajet 

10 obtenues et de I'ancienne valeur de distance affect^e au pixel en cours 
d'analyse. 

L'ordre de balayage des pixels de I'image influe sur la fiabilite des 
estimations de distance et de leurs mises ^ jour car les trajets pris en compte 
en dependent. En fait, il est soumis a une contrainte de regularity qui fait que 

15 si les pixels de I'image sont reper6s selon l'ordre lexicographique (pixels 
classes dans un ordre croissant ligne par ligne en partant du haut de Timage 
et en progressant vers le bas de Timage, et de gauche a droite au sein d'une 
ligne), et si un pixel p a 6ty analyst avant un pixel q alors un pixel p+x doit 
§tre analyses avant le pixel q+x. Les ordres lexicographique (balayage des 

20 pixels de Timage ligne par ligne de haut en bas et, au sein d'une ligne, de 
gauche d droite), lexicographique inverse (balayage des pixels de I'image 
ligne par ligne de bas en haut et, au sein d'une ligne, de droite d gauche), 
lexicographique transpose (balayage des pixels de I'image colon ne par 
colonne de gauche d droite et, au sein d'une colonne, de haut en bas), 

25 lexicographique transpose inverse (balayage des pixels par colonnes de 
droite a gauche et au sein d'une colonne de bas en haut) satisfont cette 
condition de regularity et plus g6neralement tous les balayages dans 
lesquels les lignes et colonnes, ou les diagonales sont balay6es de droite a 
gauche ou de gauche a droite. G. Borgefors preconise un double balayage 

30 des pixels de I'image, une fois dans l'ordre lexicographique et une autre dans 
l'ordre lexicographique inverse. 

La figure 2a montre, dans le cas d'une passe de balayage selon 
l'ordre lexicographique allant du coin superieur gauche au coin inferieur droit 
de i'image, les cases du masque de chanfrein de la figure 1 utilisees pour 

35 r^pertorier les trajets allant d'un pixel but plac^ sur la case centrale (case 
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index6e par 0) au pixel source en passant par un pixel du voisinage dont la 
distance a d6jd fait i'objet d'une estimation au cours du meme balayage. Ces 
cases sont au nombre de huit, disposees dans la partle sup6rieure gauche 
du masque de chanfrein. II y a done huit trajets r6pertories pour la recherche 

5 du plus court dont la longueur est prise pour estimation de la distance. 

La figure 2b montre, dans le cas d'une passe de balayage selon 
I'ordre lexicographique inverse allant du coin infSrieur droit au coin sup^rieur 
gauche de Timage, les cases du masque de chanfrein de la figure 1 utilis6es 
pour rSpertorier ies trajets allant d'un pixel but piac6 sur la case centrale 

10 (case indexee par 0) au pixel source en passant par un pixel du voisinage 
dont la distance a deja fait Tobjet d'une estimation au cours du meme 
balayage. Ces cases sont complementaires de celles de la figure 2a. Elles 
sont 6galement au nombre de huit mais dispos6es dans la partie inferieure 
droite du masque de chanfrein. II y a done encore huit trajets repertories pour 

15 la recherche du plus court dont la longueur est prise pour estimation de la 
distance. 

La transformee de distance par propagation dont le principe vient 
d'etre rappele sommairement a 6te congue d Torigine pour I'anaiyse du 
positionnement d'objets dans une image mais elle n'a pas tarde a etre 

20 appliquSe d {'estimation des distances sur une carte du relief extraite d'une 
base de donn^es d'6l6vation du terrain d maillage regulier de la surface 
terrestre. En effet, une telle carte ne dispose pas explicitement d'une 
m^trique puisqu'elle est tracee d partir des altitudes des points du maillage 
de la base de donn§es d'el^vation du terrain de la zone representee. Dans 

25 ce cadre, la transformee de distance par propagation est appliqu^e a une 
image dont les pixels sont les elements de la base de donnees d*elevation du 
terrain appartenant e la carte, c*est-a-dire. des valeurs d'altitude associees 
aux coordonnees gSographiques latitude, longitude des noeuds du maillage 
ou elles ont ete mesurees, classes, comme sur la carte, par latitude et par 

30 longitude croissantes ou decroissantes selon un tableau e deux dimensions 
de coordonnees latitude et longitude. 

Pour une navigation terrain de mobiles tels que des robots, la 
transformee de distance par propagation est utilisee pour estimer les 
distances des points de la carte d'un terrain d'evotution extraite d'une base 

35 de donnees d'eievation du terrain par rapport e la position du mobile ou une 



wo 2005/031262 



PCT/EP2004/052078 



12 

position proche. Dans ce cas, 11 est connu de tenir compte de contraintes 
statiques constituees par des zones de la carte infranchissables par le 
mobile en raison de leurs configurations accident^es. Pour ce faire, un 
marqueur de zone interdite est associ6 aux 6l6ments de la base de donn^es 

5 d*6l6vation du terrain figurant dans la carte. II signale, lorsqu*il est active, une 
zone infranchissable ou interdite et inhibe toute mise a jour autre qu'une 
initialisation, de Testimation de distance faite par la transform^e de distance 
par propagation pour r6l6ment pixel consid§r6. 

Dans le cas d'un a6ronef, les zones infranchissables evoluent en 

10 fonction du profil vertical impos6 d sa trajectoire si bien qu'une estimation de 
distance sous contraintes statiques au moyen d'une transformee de distance 
par propagation n*est pas satisfaisante. 

On propose de prendre en compte, dans la d6finition des zones 
interdites de passage, Taltitude pr6visible de Taeronef d chaque point but 

16 dont la distance est en cours d'estimation. Cette altitude previsible, qui 
depend bien evidemment du trajet emprunt6, est celle de I'a^ronef apr^s 
suivl du trajet adopts pour la mesure de distance. Uestimation de cette 
altitude prSvisible de I'a^ronef en un point but, se fait par propagation au 
cours du balayage de Timage par le masque de clianfrein d'une mani&re 

20 analogue d ('estimation de distance. Pour chaque trajet r^pertori^ allant d'un 
point but au point source en passant par un point du voisinage du point but 
dont la distance au point source et Taltitude previsible de I'aeronef ont dejd 
et^ estimSes au cours du m§me balayage, i'altitude previsible de I'aeronef 
est d^duite de la longueur du trajet et du profil vertical impose d la trajectoire 

25 de I'aeronef. Cette altitude previsible, estim6e pour chaque trajet repertoriS 
allant d'un point but dont la distance est en cours d'estimation a un point 
source plac6 a proximite de la position de I'aeronef, est utilis6e comme un 
critere de selection des trajets pris en compte dans {'estimation distance. SI 
elle est inferieure ou egale a i'altitude du point but figurant dans la base de 

30 donn^es d'elevation du terrain majoree d'une marge de s6curite, le trajet 
repertorie auquel elle est associee est ecarte et ne participe pas a la 
selection du plus court trajet. Une fois la selection du plus court trajet 
effectuee, sa longueur est prise pour distance du point but et I'altitude 
previsible pour I'adronef qui lui est associee est ^galement retenue pour 

35 {'altitude de Ta^ronef au point but. 
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La figure 3 illustre les principales Stapes du traitement effectuS 
lors de ('application du masque de chanfrein d un point but Pjj pour estimer 
sa distance pour un aeronef ayant un profil vertical de trajectoire impose. Le 
point but considers Pij est plac6 dans la case centrale du masque de 
5 chanfrein. Pour chaque point voisin Pv qui entre dans les cases du masque 
de chanfrein et dont la distance a deja 6t6 estimee au cours du meme 
balayage, le traitement consiste d : 

- lire la distance estim6e Dv du point voisin Pv (etape 30), 

- lire ('altitude Aij du point but Pg dans la base de donnSes 
10 d'6l6vation du terrain (etape 31), 

- lire le coefficient Cxy du masque de chanfrein correspondant d 
la case occup6e par le point voisin Pv(etape 32), 

- calculer la distance propagSe Dp correspondant a la somme de 
la distance estimee Dv du point voisin Pv et du coefficient Cxy 

15 affecte d la case du masque de chanfrein occupee par le point 

voisin Pv : 

(Stape 33), 

- calculer Taltitude prSvisible Ap de TaSronef aprSs 
20 franchissement de la distance D p directement d partir de la 

distance Dp si le profil vertical impost d la trajectoire de 
I'aeronef est deftni en fonction de la distance parcourue PV(Dp) 
et prend implicitement en compte le temps de parcours ou 
indirectement par TintermSdiaire du temps de parcours si le 
25 profil vertical impose a la trajectoire de TaSronef est defini par 

une Vitesse de changement d'altitude (etape 34), 

- comparer I'altitude previsible Ap obtenue avec celle Aj.j du point 
but Pij tiree de la base de donnees d'elevation du terrain 
augmentee d'une marge de s6curit6 A (etape 35), 

30 - eliminer la distance propagSe Dp si I'altitude previsible Ap est 

inferieure ou egale a celle Ay du point but Pij tiree de la base 
de donnSes d'Stdvation du terrain et majoree par la marge de 
sScurite A (etape 36), 

- si I'altitude previsible Ap est supSrieure d celle Aij du point but 
35 Pij majorSe par la marge de security A, lire la distance Dy deja 
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affect6e au point but consid6r§ Pij ( 6tape 37) et la comparer d 

la distance propagde Dp (6tape 38), 
- 6liminer la distance propag6e Dp si elle est sup6rieure ou 

Sgale d la distance Dy ddja affect6e au point but considere Pij, 
5 et 

. remptacer la distance Dy dej^ affect6e au point but consider6 

Pg, par la distance propag^e Dp si cette derniere est inferieure 

(etape 39). 

La figure 4 illustre les principales 6tapes d'une variante du 

10 traitement effectue lors de I'application du masque de chanfrein a un point 
but Pij pour estimer sa distance pour un aeronef ayant un profil vertical de 
trajectoire impose. 

Cette variante diff6re dans la maniere d'6laborer {'altitude 
previsible Ap de Taeronef. Elle suppose que {'altitude previsible pour {'a6ronef 

15 en chaque point de la base de donn6es d*§l6vation du terrain calculee en 
fonction du profil vertical impos6 d sa trajectoire et ^ partir de la longueur du 
trajet selectionn6 pour la mesure de distance n*est pas consid§r6e comme 
une variable passagdre, ce que permet le traitement d^crit relativement d la 
figure 3, mais est m§morisee, au meme titre que Testimation de distance. 

20 Dans cette variante, les dtapes 30, 31 de lectures de la distance estim^e Dv 
du point voisin Pv et de Taltitude Ay du point but Py dans la base de donn6es 
d'6l6vation du terrain sont completees par une etape 40 de lecture de 
{'altitude previsibie Apv pour {'aeronef au point voisin Pv, et ie calcul de 
{'altitude previsible Ap se fait (etape 34') par sommation de {'altitude previsible 

25 Apv au point voisin Pv et de la variation d'altitude sur la distance separant (e 
point voisin Pv du point but Pg due au profil verticai impose d ia trajectoire de 
{'aeronef. 

La memorisation des aititudes previsibies pour {'a6ronef {orsqu'iF 
atteint ies differents points de la carte qui iui sont accessibles permet 

30 d'etablir, en {eur soustrayant les altitudes des points de la carte tirees de ia 
base de donn^es d'6(6vation du ten^in, une carte des possibilites maximaies 
de hauteurs de survol de Taeronef representant les hearts previsibies par 
rapport au terrain en strates de couleurs. Une telle carte facilite pour 
('equipage de I'a^ronef le choix d'une trajectoire r^aliste pr^sentant la 

35 meilleure garde au sol. 
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Comme indique pr6cedemment, Testimation des distances des 
differents points de la carte se fait en appliquant un traitement par masque 
de chanfrein tel que ceux qui viennent d'etre d6crits relativement aux figures 
3 et 4, ^ Tensemble des pixels de I'image formee par ies 6l6ments de la base 
5 de donnSes d'^l^vation du terrain appartenant d la carte, pris 
successivement selon un balayage r^gulier comportant un minimum de deux 
passes r6alis6es dans des ordres inverses. 

La figure 5 illustre Ies principales etapes d'un exemple de 
processus global permettant Testimation des distances de {'ensemble des 
10 points d'une carte de relief pour un mobile soumis ^ des contraintes 
dynamiques. 

La premiere ^tape 50 du processus est une initialisation des 
distances affectees aux differents points de la carte consideres comme Ies 
pixels d'une image. Cette initialisation des distances consiste, comme 

16 indiqu6 precedemment, a attribuer une valeur de distance infinie, a tout le 
moins sup6rieure a la plus grande distance mesurable sur la carte, pour tous 
Ies points de la carte consideres comme des points but, a Texception d'un 
seul consider^ comme la source de toutes Ies distances auquel est attribu6 
une valeur de distance nulle. Ce point source est choisi S proximity de la 

20 position instantan^e du mobile sur la carte. 

Les 6tapes suivantes 51 d 54 sont des passes d'un balayage 
regulier au cours desquelles le masque de chanfrein est appliqu6 
successivement et ^ plusieurs reprises d tous les points de la carte 
consid6r6s comme les pixels d'une image, I'application du masque de 

25 chanfrein d un point de la carte donnant une estimation de la distance de ce 
point par rapport au point source, par execution d'un des traitements d^crits 
relativement a la figure 3 ou d la figure 4. 

La premiere passe de balayage (etape 51) se fait dans Tordre 
lexicographique, les pixels de Timage etant analyses ligne par ligne du haut 

30 vers le bas de Timage et de gauche S droite au sein d'une meme ligne. La 
deuxieme passe de balayage (§tape 52) se fait dans I'ordre lexicographique 
inverse, les pixels de I'image 6tant toujours analyses ligne par ligne mais du 
bas vers le haut de I'image et de droite a gauche au sein d'une ligne. La 
troisieme passe de balayage (etape 53) se fait dans I'ordre lexicographique 

35 transpose, les pixels de {'image etant analyses colonne par colonne de la 
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gauche vers la droite de Timage et de haut en bas au sein d'une meme 
colonne. La quatridme passe de balayage (§tape 54) se fait dans Tordre 
lexicographique transpose inverse, les pixels de Timage etant analyses 
colonne par colonne mais de la droite vers la gauche de rimage et de bas en 
5 haut au sein d'une meme colonne. 

Ces quatre passes ( 6tapes 51 a 54) sont rep^tees tant que 
rimage de distance obtenue change. Pour ce faire, le contenu de Timage de 
distance obtenu est m6moris6 (etape 56) apres chaque s§rie de quatre 
passes ( etapes 51 a 54) et compare avec le contenu de Timage de distance 

10 obtenu a la serle precedente (etape 55), la boucle n'6tant brisee que lorsque 
la comparaison montre que le contenu de rimage distance ne varie plus. 

En theorie, deux passes de balayage selon I'ordre lexicographique 
et I'ordre lexicographique Inverse peuvent suffire. Dependant la pr6sence de 
zones interdites de passage de forme concave peut provoquer, dans le 

15 phenom6ne de propagation des distances, des angles morts renfermant des 
pixels pour lequel I'appiication du masque de chanfrein ne donne pas 
d'estimatlon distance. Pour diminuer ce risque d'angle mort, il y a lieu de faire 
varier la direction du phenomene de propagation distance en faisant varier la 
direction du balayage d'ou le doublement des passes avec une transposition 

20 des ordres de balayage correspondant d une rotation de l*image de 90"". Pour 
une encore meilleure Elimination des angles morts, on peut procSder d des 
series de huit passes : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
rimage. chaque ligne etant parcourue de gauche d droite, 

25 - une deuxieme passe effectu6e ligne par ligne de bas en haut 

de rimage, chaque ligne 6tant parcourue de droite d gauche, 

- une troisieme passe effectu6e colonne par colonne de gauche 
d droite de rimage, chaque colonne etant parcourue de haut 
en bas, 

30 - une quatri§me passe effectuee colonne par colonne de droite ^ 

gauche de I'image, chaque colonne etant parcourue de bas en 
haut, 

- une cinqui^me passe effectuee ligne par ligne de haut en bas 
de rimage, chaque ligne Etant parcourue de droite d gauche, 
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- une sixi6me passe effectu6e ligne par ligne de bas en haut de 
rimage, chaque ligne 6tant parcoume de gauche a droite, 

- une septieme passe effectuee colonne par colonne de droite ^ 
gauche de I'image, chaque colonne 6tant parcoume de haut en 

5 bas, et 

- une huiti^me passe effectuee colonne par colonne de gauche 
gi droite de Timage, chaque colonne Stent parcourue de bas en 
haut. 

II est possible d'Introduire dans la serie de passes de balayage 
10 d'autres types de passes de balayage d6duits des passes precedentes en 
faisant jouer par les diagonales de rimage, les roles jou6s pr6cedemment 
par les lignes et colonnes de rimage. Cela revient S appliquer les passes de 
balayage decrites prec6demment S une image tournee de 45°. De mani6re 
generate, plus les passes d'une s§rie sont vari6es plus le risque d'angle mort 
15 diminue. 
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REVENDICATIONS 

1, Precede d'estimation des distances des points d'une carte 

5 extraite d'une base de donnees d'6l6vation du terrain, pour un mobile soumis 
a des contraintes dynamiques lui interdisant certaines zones de la carte dites 
zones interdites de passage dont la configuration varie en fonction du temps 
de parcours du mobile ; la base de donnees d*elevation du terrain renfermant 
un ensemble de points repdrds par une altitude, une latitude et une longitude 

10 maillant le terrain d'evolution du mobile ; ledit proc^dd mettant en oeuvre une 
transformee de distance operant par propagation sur Timage constituee par 
les elements de la base de donnees d'elevation du terrain correspondant ^ la 
carte, disposes en lignes et colonnes par ordres de valeurs de longitude et 
de latitude ; la transform§e de distance estimant les distances des diffSrents 

15 points de I'image par rapport d un point source place A proximite du mobile, 
en appliquant, par balayage, un masque de chanfrein aux differents points de 
rimage ; i'estimation de distance d'un point, par application du masque de 
chanfrein d ce point dit point but s'effectuant en r6pertoriant les differents 
trajets allant du point but au point source et passant par des points du 

20 voisinage du point but qui sont couverts par le masque de chanfrein et dont 
les distances au point source ont ete prealablement estim^es au cours du 
meme balayage, en determinant les longueurs des differents trajets 
repertories par sommation de la distance affectee au point de passage du 
voisinage et de sa distance au point but extraite du masque de chanfrein, en 

25 recherchant le trajet le plus court parmi les trajets rdpertorids et en adoptant 
sa longueur comme estimation de la distance du point but ; une valeur de 
distance superieure d la plus grande distance mesurable sur I'image 6tant 
initialement attribuee, en d§but de balayage, ^ tous les points de I'image sauf 
au point source, origine des mesures de distance, auquel est affectSe une 

30 valeur de distance nulle ; ledit proc^dS etant caract^risS en ce que les 
longueurs des trajets repertories, lors de i'application du masque de 
chanfrein a un point but, en vue de la recherche du trajet le plus court, sont 
traduites en temps de parcours pour le mobile et en ce que les trajets 
repertories dont les temps de parcours pour le mobile sont tels que le point 

35 but appartiendrait ^ une zone interdite de passage au moment ou le mobile 
Tatteindrait, sont exclus de la recherche du trajet le plus court. 
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2. Procede selon la revendication 1, applique d un a^ronef ayant 
un profil vertical de vol ^ respecter determinant revolution de son altitude 
instantan^e, caracterlse en ce que, aux longueurs des trajets repertoriees 

5 lors de Vapplication du masque de chanfrein a un point but, sont associ^es 
les vaieurs previsibles des altitudes instantan§es qu*aurait Taeronef en 
atteignant le point but par ces trajets tout en respectant le profil vertical de 
vol impose, et en ce que les trajets repertories associes a des vaieurs 
previsibles d'altitude inf^rieures ou egales a celle du point but donn^e par la 
10 base de donnees d'el^vation du terrain et augment^e d'une marge de 
protection sont exclus de la recherche du trajet le plus court,. 

3. Proc6d6 selon la revendication 2, caracteris6 en ce que 
Testimation de distance operee par propagation sur l*image constitute des 

15 elements de la base de donnees d'6l6vation du terrain correspondant a la 
carte est doublte d'une estimation de Taltitude previsible pour Tatronef au 
droit des difftrents points de I'image en supposant qu*il suive le trajet le plus 
court seiectionne pour Testimation de distance et qu'il respecte le profil 
vertical de vol impost. 

20 

4. Proctdt selon la revendication 3, caracterist en ce que les 
altitudes des difftrents points de la carte sont soustraites des estimations 
des altitudes previsibles pour Tatronef en ces points pour obtenir des ecarts 
par rapport au sol. 

25 

5. Procede selon la revendication 4, caracterise en ce que les 
ecarts par rapport au sol sont affiches sur la carte en strates de couleurs. 

6. Precede selon la revendication 1, caracterise en ce que la 
30 transformee de distance par propagation balaie les pixels de Timage 

constituee des elements de la base de donnees d'eievation du terrain 
correspondant a la carte, en plusieurs passes successives selon des ordres 
differents. 
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7. Procede selon la revendication 6. caracteris^ en ce que la 
transformSe de distance par propagation balaie les pixels de rimage 
constitute des fitments de la base de donntes d'tltvatlon du terrain 
appartenant a la carte, en plusleurs passes successives selon des ordres 

5 difftrents et de manitre rtpetde jusqu'a ce que les estimations de distance 
obtenues se stabllisent. 

8. Proc6de selon la revendication 6, caract6rise en ce que la 
transformee de distance par propagation balaie les pixels de rimage 

10 constituee des el6ments de la base de donn6es d'6l6vation du terrain 
correspondent a la carte, en plusieurs passes successives selon des ordres 
differents dont I'ordre lexicographique, I'ordre lexicographique inverse, Tordre 
lexicographique transpose et I'ordre lexicographique transpose inverse. 

15 9. Proc6de selon la revendication 6, caracteris6 en ce que la 

transformee de distance par propagation balaie les pixels de I'image 
constitute des Elements de la base de donntes d'tltvation du terrain 
correspondant a la carte, en une strie de quatre passes repttte jusqu'a 
stabilisation des estimations de distances : 

20 - une premiere passe effectuee ligne par llgne de haut en bas de 

rimage, chaque ligne ttant parcourue de gauche d droite, 

- une deuxieme passe effectute ligne par ligne de bas en haut 
de rimage, chaque ligne ttant parcourue de droite a gauche, 

- une troisieme passe effectute colonne par colonne de gauche 
25 d droite de I'image, chaque colonne etant parcourue de haut 

en bas, et 

- une quatrieme passe effectuee colonne par colonne de droite ci 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut. 

30 

10. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constituee des elements de la base de donnees d*tltvation du terrain 
correspondant a la carte, en une serie de huit passes repetee jusqu'a 
35 stabilisation des estimations de distances : 
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- une premiere passe effectuSe ligne par ligne de haut en bas de 
rimage, chaque ligne Stant parcourue de gauche d droite, 

- une deuxidme passe effectude ligne par ligne de bas en haut 
de rimage, chaque ligne 6tant parcourue de droite d gauche, 

5 - une troisiSme passe efYectu6e colonne par colonne de gauche 

a droite de Timage, chaque colonne 6tant parcourue de haut 
en bas, 

- une quatridme passe effectu6e colonne par colonne de droite d 
gauche de I'image, chaque colonne 6tant parcourue de bas en 

10 haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de rimage, chaque ligne etant parcourue de droite a gauche, 
une sixieme passe effectuee ligne par ligne de bas en haut de 
rimage, chaque ligne 6tant parcourue de gauche ^ droite, 

15 - une septi6me passe effectuee colonne par colonne de droite d 

gauche de Timage, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe efFectu6e colonne par colonne de gauche 
d droite de Timage, chaque colonne etant parcourue de bas en 

20 haut. 

1 1 . Proc§d6 selon la revendication 6, caract^risd en ce que la 
transformee de distance par propagation balaie les pixels de Timage 
constituee des elements de la base de donn^es d'el§vation du terrain 
25 appartenant a la carte, en plusieurs passes successives selon des ordres 
differents dont certains consistent en un balayage de I'image par diagonales, 
d*un bord a Tautre et, au sein d'une diagonale, d*une extremite a I'autre. 
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